Spin-liquid model of the sharp resistivity drop in Lai^s^Bao^u^CuOi. 
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We use the phenomenological model proposed in our previous paper [Phys. Rev. Lett. 98, 237001 
(2007)] to analyse the magnetic field dependence of the onset temperature for two-dimensional 
fluctuating superconductivity T**{H). We demonstrate that the slope of T**{H) progressively goes 
down as H increases, such that the upper critical field progressively increases as T decreases. The 
quantitative agreement with the recent measurements of T**{H) in Lai,ssBao.i25Cu04 is achieved 
for the same parameter value as was derived in our previous publication from the analysis of the 
electron self energy. 
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Recent experiments on Lai^xBaxCuOi at x = 1/8 
[l[ revealed a complex hierarchy of energy scales in this 
material. It displays a charge ordering transition at 
Tco = hAK , a spin ordering transition at Tspin — 42K 
with a subsequent one order of magnitude drop in the 
in-plane resistivity, the Berezinskii-Kosterlitz-Thouless 
(BKT) transition to a two-dimensional superconductiv- 
ity at Tbkt — lOif, a crossover from 2D to 3D regime 
around \QK , and a transition to a true 3D superconduc- 
tivity at AK. This hierarchy is summarized and discussed 
in detail in f3|. 

It turns out that the temperature T** where the resis- 
tivity crossover occurs is sensitive to the c-axis magnetic 
field which separates this phenomenon separately from 
the spin ordering. In this paper, we address the issue of 
this crossover. The measurements performed in a mag- 
netic field [H revealed that (i) T** marks the onset of 
fluctuational diamagnetism, and (ii) T** decreases with 
the field. These two effects and the fact that the resistiv- 
ity sharply drops T** are consistent with the idea that 
T** marks the onset of a fluctuational pairing regime 
without (quasi-) long-range superconducting order. The 
details of the system behavior near T**, however, depend 
on the underlying model. The authors of 0] consid- 
ered a model of weakly coupled parallel superconducting 
stripes. Within this model, T** is the temperature at 
which the inter-stripe coupling becomes strong, and a 
vortex liquid is formed. 

We propose another explanation, based on the model 
with a flat Fermi surface in the antinodal regions near 
(0, tt) and (tt, 0) points in the Brillouin zone Q. Fermions 
in these regions form two quasi-lD spin liquids coupled 
by Josephson-type interaction. In this model, the pair- 
ing amplitudes in the antinodal regions are developed at 
T* ^ T** due to the attractive interactions in the spin- 
liquid state, however, phase fluctuations at T >> T** 
are effectively one-dimensional, and are pinned by the 
defects. At T**, the Josephson coupling becomes suffi- 
ciently strong to lock the relative phase of the two or- 
der parameters at tt, and the system response becomes 
two-dimensional. This leads to depinning of the phase 




FIG. 1: The plot of T**{H). The points are the data from 
Ref. [ll], the solid curve is the exponential fit by Eq. ([l]). 



fluctuations resulting in the drop in the resistivity. Still, 
because of vortices in the 2D regime, the (quasi)-long- 
range superconducting order develops only at a smaller 
Tc < T**. 

Just like the model of parallel stripes our model of 
"crossed stripes" near (0,7r) and (tt, 0) explains qualita- 
tively the resistivity drop, the absence of fluctuational 
diamagnetism above T**, and the sensitivity of T** to 
a magnetic field. [H]- However, the measurements of 
T**{H) put an additional constraint on the theory - not 
only T** decreases with the field, but \dT** /dH\ also de- 
creases as H goes up, i.e., at very low T, the critical field 
below which the system response is two-dimensional, be- 
comes very large. The data for H < 9T can be well fitted 
by the exponential dependence (see Fig. 1): 



T**{H) 

r**(o) 
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For such T**{H), \dT**/dH\ exponentially decreases as 
H increases. If this trend continued to higher H, the crit- 
ical field Hc2{T) defined as T**{Hc2) = T would become 
infinite at T = 0. 

The H dependence of T** for Josephson-coupled 
stripes running parallel to each other in the 2D plane. 
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i.e., for the same model as in Ref. Q was considered by 
Carr and one of us It was found that the sfope of 
dT** /dH increases with decreasing T, and ifc2 remains 
finite at T = 0, in quahtative disagreement with the 
data. We demonstrate below that our model of crossed 
stripes located near (0, tt) and (tt, 0) yields the behavior 
of T**{H) in a good agreement with the measurements. 
Thus we show that the slope of dT** /dH decreases with 
increasing H for any value of the scaling dimension d of 
the superconducting order parameter. To achive a quan- 
titative agreement with the experimental fit ([1]) we have 
to set d « 1/2. We have to remind the reader that in [3] 



the same value of d was postulated on the basis of anal- 
ysis of the electron self energy. This gives an important 
check for self-consistency of the theory. 

We associate T**{H) with the instability of a 2D pair- 
ing susceptibility in the random phase approximation 
(RPA). Fluctuations beyond RPA transform the instabil- 
ity into a crossover 0] . In zero field, the RPA expression 
for the susceptibility reads, in momentum space 

) Xo{kx) Xo{ky) (2) 
where Xo(fc) is the ID static pairing susceptibility 0: 



sin^d T^{l-d) (^^^ 



-2+2d 



T{d/2 + ivq/4TTT) 



r(l - d/2 + ivq/47rT) 
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Here r(...) are F— functions, d < 1 is the scaling dimen- 
sion of the superconducting order parameter, v is the 
velocity of the phase mode, and A is the ultraviolet cut- 
off. The last term in xo can be neglected as we will 
only consider T <C A, when the first term in ([3]) dom- 
inates. Parameters v and d are free parameters of our 
theory and should be extracted from the experiments in 
the T region where the superconducting phase fluctua- 
tions are essentially one-dimensional (that is, at T below 
the spin gap, but larger than T**). In \^ we found that 
the best agreement with the photoemission experiments 
is obtained when d « 1/2. As we will see, this value is 
also favored by the observed T**{H) dependence. 

Taking a Fourier transform over kx, but leaving ky 
intact, we obtain from 

Xky{x-xi)= (4) 
Xo(a; - xi) / dx'xQ{ky)xo{x - x')xky{x' - xi) 



In a magnetic field, ky 



kii 



Hx' (we set 2e/c = 1) 



Setting ky = and xi = 0, we obtain integral equation 
for xix) = Xky=o{x) in the form 

X{x) = xo{x) + I dx'xoix - x')x{x')xoiHx') (5) 



where xq{Hx') is given by ([3]) for fc = Hx', and Xo{x) 
is the Fourier transform of Xo(^)- The temperature 
T**{H) is the one at which x{x) diverges. 

Weak fields. Consider first the case when the mag- 
netic field is weak, i.e., T**{H) = T**(0)(1 - 5T), and 
5T ^ 1. A simple analysis shows that the parametrical 
condition for a weak field is v^H/T « 1. Expanding 



Xo{Hx') in H, we obtain from Q 
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V'^^^ (d/2) - V^^^ {l-d/2) 
F2 (d/2) 



(7) 



and V'^-^H^) is the derivative of the diGamma function. 

Substituting ^ into ([5]), we obtain an integral equa- 
tion for x{x) in the form 

X{x) = Xo{x) + J dx'xoix - a;')x(a;')Xo(0) 
v^H^ 



-rxo{o)Ad 



dx'xoix -x')x{x')(x'f (8) 



where xo(0) ~ Xoik — 0) . Taking Fourier transform back 
to momentum space (x — > fcj; — k), and integrating by 
parts, we re- write the integral equation for x as 

X(fc) [1 - J\o{k)xom - J\o{k)xo{0) 



- xoik) 
This can be re-expressed as 

d 



(9) 



e + cik - C2^j x(fc) = Xo(fc) 
1 - (T**(0)/T)4-4<i^ 



(10) 



where e = 1 - (T**(0)/T)4-4^ ci = Adv'^/{nT)^, 

C2 



AdV^H^/inTy, and we defined r**(0) 
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(A/27r) (SdJ)i/(2-2d). This agrees with the zero-field 
transition temperature in [3|. Expanding now in the 
eigenvalues of the difl^erential equation as 



where Xn{k) are the solutions of 

Clk^ -C2-Q^\Xn{k)=€nXn{k), (12) 



we obtain 



,(0) 



(13) 



The eigenvalues of Eq. (fT2|l can be easily obtained as 
(fT^ can be re-expressed as a harmonic oscillator 



1 d'^Xn{k) Muj'^k"' 



2M 9fc2 



X«(fc) =e«X„(fc) (14) 



where uj^ = 4ciC2 and M^^ = 2Ad('y/7rr)2. The eigen- 
functions of (|14p are e„ = ^^(n + 1/2), the lowest one 
is eo = t^/2 = Av'^H/{TrT)'^. From ([131), the insta- 
bility in the field occurs when e -I- eo — 0, i.e, when 
T = T**{H) = T**(0)(1 - ST), where 



ST 



1 



4(1 -d) (7rT**(0))2 



(15) 



We see that at small fields, T**{H) decreases linearly 
with H . The linear dependence at small fields is also 
present in the model of parallel stripes Q • If we formally 
extrapolate the small-field result to T = 0, we obtain the 
upper critical field 



Hir{T = 0) = [-\ (Ji3,)V(i-<i)L^ 
V I Ad 



(16) 



The actual Hc2{T = 0) is somewhat smaller in the model 
of parallel stripes |J], but, as we will see, is much larger 
than in our model of crossed stripes. 

Strong fields. Consider now the opposite limit of 
vanishing T, when v^H/T » 1, i.e., the expansion in 
the field is no longer possible. In this limit, we have from 




Xo{Hx') = 



Hx 



f\2-2d 



(17) 



where 



Bd = (8/A)2sin(7rd) T'^{1 - d){v'^ /AA^ 

= Bd{2A/vf-^'^ (r2(l - d/2)/r2(d/2)). (18) 

Instead of Eq. ([9]) , we now have 

J(k-q)x' 



X{k) = Xo(fc) 



1 + J 



2 B2d-1 
Jj2~2d 



dqxiq) / dx' 



„/|2-2d 



(19) 



Using 



,e*(fe-9)^' _ r(2(i- l)sin7rd 



and introducing 



\k\' 



and d = (1 + e)/2, we obtain from ([T 



(21) 



m = 1+ 



J2b2 cos(7re/2) r(e) 



H 



q\^-'\k-q\' 



(22) 



It is convenient to re-express this equation in the operator 
form, as Lx{k) = 1, and expand in the eigenfunctions of 
the operator L, which we label as Xm{k)- We get 



(0) 



X{k) = V Xm(fc) 
— i — Am. 



(23) 



where am"' are constants. The eigenvalues A™ are the 
solutions of 



LXm{k) = (1 - \n)Xm{k) 



(24) 



where 



LXm{k) = Xm{k)- 



j'^Bj cos7re/2r(e) 



Xm{q) 



\q\'-'\k-q\' 
(25) 

Eq. ((25)) was studied in the context of non-BCS super- 
conductivity (with frequency instead of momentum) 
A similar equation has been studied in the content of 
superconductivity in graphene Q. For e > 0, the nor- 
malized solution of with the largest eigenvalue is 



Xmik) 



(26) 



and the eigenvalue is 



J'B 



1 



(1 - e/2)(sin7re/4)2 



(27) 



The critical field Hc2 {T = 0) is determined from Ao = 1 
and is given by 



Hc2{T = 0) = 

In explicit form, we have 



2o2,T, 1 1/(1-0 



(28) 
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Hc2{T = 0) = {JBd) 



i/(i-rf) f ^ 
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FIG. 2: The ratio 77^2 (T = 0)/H, 
scaling dimension d, from Eg. ([29 
rithmicaly at d ^ 0.5. 



extr 
c2 



(T = 0) as a function of 
The ratio diverges loga- 



One can easily make sure that the actual idc2 [T = 0) is 
much larger than iJ^f*''(T = 0) for all d < 1/2 for which 
our computational scheme is applicable. Furthermore, as 
d approaches 1/2, Hc2{T = 0) tends to infinity because 
^'(e) diverges at vanishing e = 2d — 1 as w 8/e^. The 
plot of the ratio Hc2{T = 0)/H^^*''{T = 0) is presented 
in Fig. 2. 

For d < 1/2, the analysis has to be modified to ac- 
count for the divergence at g = in the r.h.s. of 
The expected result is that Hc2 becomes infinite at zero 
temperature. The divergence is power-law for e < 0, and 
logarithmical at e = 0. In the latter case. 



Xo{Hx') 



Hx' 



(30) 



and the RPA equation for x(x) in the real space becomes 

x(x) = -21ogTN-^^^^§^ J 0jx(x') \og{T\x-x'\) 

(31) 

With the logarithmic accuracy, we can approximate 
log (|a; - x'\) « 0ix - x') log \x\ + e{x' - x) log \x'\ (32) 



Substituting into ([31]) . we re- write it as a differential 
equation 



d'a + X = -2d', log {T\e< - 1|) (33) 



where C — ^og\x\. The analysis of this equation shows 
that the susceptibility diverges at if = Hc2{T) oc \logT\. 
This is equivalent to T**{H) oc exp—H/Ho, in agreement 
with Eq. (|T]). We see therefore that the high field depen- 
dence is well captured by our model with d w 1/2 - the 
same as we used in the previous work [3| to fit the normal 
state self-energy. 

To summarize, we analyzed the behavior of T**(H) 
(or, equivalently Hc2(T)) in the model of two one- 
dimensional spin liquids near (0, tt) and (tt, 0) coupled 
by Josephson-type interaction. For weak fields we found 
that T** decreases linearly with H. Extrapolating this 
dependence down to zero temperature yields the ex- 
trapolated field H'^2*^{T = 0). Considering the strong 
fields we found that the actual Hc2{T — 0) is always 
larger than the extrapolated value. The ratio Hc2{T = 
0)/H^2*^{T — 0), characterizing the convexity of the 
_ffc2(r)-curve, increases when d decreases and becomes 
infinite at d < 1/2. This convex behavior is consistent 
with the data, and has to be contrasted with the concave 
behavior for the model of parallel stripes. As a further 
evidence in support of our model, we found that the ex- 
perimental Hc2{T) are well described by the theoretical 
formula with the scaling dimension of the ID supercon- 
ducting order parameter d w 1/2. The same d provides 
the best fit to the photoemission data, as we argued ear- 
lier (3] . We think that all these give our model a consid- 
erable advantage in treating Lai,s5Bao.i25CuOi. 
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